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ABSTRACT 

We present and employ a new kinematical approach to cosmological 'dark energy' 
studies. We construct models in terms of the dimensionless second and third derivatives 
of the scale factor a{t) with respect to cosmic time t, namely the present-day value of 
the deceleration parameter qo and the cosmic jerk parameter, j(t). An elegant feature 
of this parameterization is that all ACDM models have j{t) = 1 (constant), which 
facilitates simple tests for departures from the ACDM paradigm. Applying our model 
to the three best available sets of redshift-independent distance measurements, from 
type la supernovae and X-ray cluster gas mass fraction measurements, we obtain clear 
statistical evidence for a late time transition from a decelerating to an accelerating 
phase. For a flat model with constant jerk, j(t) = j, we measure qo = —0.81 ± 
0.14 and j = 2.16l^;^5, results that are consistent with ACDM at about the la 
confidence level. A standard 'dynamical' analysis of the same data, employing the 
Friedmann equations and modeling the dark energy as a fluid with an equation of 
state parameter, w (constant), gives 0^ = 0.306to;o4o ^^'^ ^ = ^l-15lo;i8, also 
consistent with ACDM at about the Icr level. In comparison to dynamical analyses, 
the kinematical approach uses a different model set and employs a minimum of prior 
information, being independent of any particular gravity theory. The results obtained 
with this new approach therefore provide important additional information and we 
argue that both kinematical and dynamical techniques should be employed in future 
dark energy studies, where possible. Our results provide further interesting support 
for the concordance ACDM paradigm. 

Key words: cosmology: observations - cosmology:cosmological parameters - cosmol- 
ogyitheory - supernovae - x-ray clusters - dark energy 



1 INTRODUCTION 

The field of cosmology has made unprecedented progress 
during the past decade. This has largely been driven 
by new observations, including precise measurements of 
the spectru m of cosmic microwave background (C MB) 
anisotropies (jSpergel et al.l l2003l : ISpergel et all |2006| . and 
references there in), the dis t ance- r edshift relation to type 
la supernov ae (iRiess et all Il998l: iPerlmutter et"al] Il999l : 



Eisenstein et al.ll2005h. Lvman-a forest studies (ICroft et al.l 


20021: Viel et al.ll2004 ISeliak et al.l 


20051) and weak lensing 


surveys (Ivan Waerbeke et al. 200C 


: iHoekstra et al.l |2002|: 



measure- 



Astier et al 



gala xy clusters 



iKnop et al.1 l2003l : iRiess et al.l 12004 ? 

the distance - redsh if t rela t ion to X-ra y 
lAUen et al.l |2002| . 12004 lEttori et all |2003| ). measure- 
ments of the mean matter density and amplitude of 
matt er fluctuation s from X-ray clusters (iBorgani et al.l 
20011: lAUen et al.1 [200^: iReiprich fc BohringeJ l200a 



Schuecker et all l2003l : IVoevodkin fc VikhlininI |2004| ). mea 
surements of the matter power spec tr um from galax; 
redshift surveys (|Tegmark et al.1 |2004| : ICole et al.1 1200 
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ments of the Integrated Sachs-Wolfe effect (|Scranton et al.l 
I2OO3I : iFosalba et atlbOOSl ) . 

These and other experiments have lead to the defini- 
tion of the so-called concordance ACDM cosmology. In this 
model, the Universe is geometrically flat with only ~4 per 
cent of the current mass-energy budget consisting of normal 
baryonic matter. Approximately 23 per cent is cold dark 
matter, which interacts only weakly with normal baryonic 
matter but which clusters under the action of gravity. The 
remaining ~ 73 per cent consists of smoothly distributed 
quantum vacuum energy (the cosmological constant), which 
pushes the Universe apart. This combination of matter and 
vacuum energy leads to the expectation that the Universe 
should undergo a late time transition from a decelerating to 
an accelerating phase of expansion. Late-ti me acceleration of 
the Universe is now an observed fa ct (e.g. iRiess et al. I I2OO4I : 
I Allen et al]|2004l : lAstier et al.ll2005h . A transition from a de- 
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celerating phase to a late-time accelerating phase is required 
to explain both these late-time acceleration measurements 
and the observed growth of structure. 

Despite the observational success of the concordance 
ACDM model, significant fine tuning problems exist. In 
particular, difficulties arise in adjusting the density of 
the vacuum energy to be a non-zero but tiny number, 
when compared with the value predicted by standard the- 
oretical calculations, and with explaining why the cur- 
rent matter and vacuum energy densities are so simi- 
lar (the 'cosmic coincidence' problem). For these reasons, 
amongst others, a large number of alternative cosmologi- 
cal models have been proposed. These include models that 
introduce new energy components to the Universe - so 



JCaldwell et al. 19981: 


Zlatev et al. 19991: Cooeland et al.l 


19981: ISteinhardt et al 


Il999l: iBarreiro et al.l l2000t). K- 


essence ( Chiba et al.l 2000l: lArmendariz-Picon et al.l |2000|. 


200ll). tachvon fields (Bagla et al.l 


2003; Cooeland et al.1 


2005|~) and Chaplvein eas models (Kamenshchik et al.ll200ll: 



Bento et al.ll2002 ) . Other possibilities include modified grav- 



ity th eories, motivated by e.g. the existence of extra dimen- 



2OO2I: Maartens & Maierotto. ,200a: .Guo ct al. 


.2006,1 or 


other modifications of General Relativitv (ICaDOzziello et al. 


2002 


: ICarroU et al.1 12004 


; IVoUickl I2OO3I: Carroll et al. 


2005 


: Mena et all I2OO6I: 


Navarro & Van Acoleven 2003: 


Noiiri & Odintsovl I2OO6I) . which can also lead to late-time 



cosmic acceleration. The simplicity of the concordance 
ACDM model, however, makes it highly attractive. A central 
goal of modern observational cosmology is to test whether 
this model continues to provide an adequate description of 
rapidly improving data. 

Most current analyses of cosmological data assume Gen- 
eral Relativity and employ the mean matter density of the 
Universe, fim, and the dark energy equation of state w as 
model parameters. Such analyses are often referred to as 
'dynamical studies', employing as they do the Friedmann 
equations. Other dynamical analyses employ modified Fried- 
mann equations for a particular gravity model. However, a 
purely kinematical approach is also possible that does not 
assume any particular gravity theory. Kinematical models 
provide important, complementary information when seek- 
ing to understand the origin of the observed late-time accel- 
erated expansion. 

In a pioneering study, iRiess et al.l (|2004l ') measured a 
transition from a decelerating to accelerating phase using a 
simple linear parameterization of the deceleration param- 
eter q{z), where q{z) is the dimensionless second deriva- 
tive of th e scale factor, a(t) , with r espect to cosmi c time . 
Rec entlv. IShapiro fc Turner (l2005h . ICong fc Wa"n3 (|2006l ) 
and lElgarov &: Multamakil ( 2006t ) have employed a variety 
of other parameterizations, constructed in terms of q{z), to 
study this transition. However, since the underlying physics 
of the transition are unknown, the choice of a particular pa- 
ramet erization for q{z) is quite arbitrary. IShapiro fc Turnen 
l|2005h applied a pr incipal component analysis of q{z) to the 
supernovae data of iRiess et al.l (|2004l ') and found strong ev- 
idence for recent, changing acceleration but weak evidence 
for a decelerated phase in the past (i .e. weak evidence for a 
transi tion between the two phases). lElgarov fc Multamaki 
l|2006l ) employed a Bayesian analysis to the Riess et al 



(|2004|) data and the more recent SNLS supernovae sample 
of lAstier etHI (|2005l '). obtaining a similar result. 

In this paper we develop an improved method for study- 
ing the kinematical history of the Universe. Instead of using 
^ aramctc rizations con structed in terms of q{z) , we follow 
Blandford et al.l (|2004l ) and introduce the cosmic jerk, j{a), 
the dimensionless third derivative of the scale factor with 
respect to cosmic time. (Here a is the cosmic scale factor, 
with a — 1/1 + z.) The use of the cosmic jerk formalism 
provides a more natural parameter space for kinematical 
studies. Our results are presented in terms of current decel- 
eration parameter qo and j{a), where the latter can be ei- 
ther constant or evolving. We apply our method to the three 
best current kinematical data sets; the'gold' sa mple of type 
la sup ernovae (hereafter SNIa) measurements of lRiess et al.l 
(|2004 ). the SNIa data from the first year of the S upernova 
Legacy Survey (SNLS) project (jAstier et al.ll2005l'l. and the 
X-ray galaxy cluster distance measurements of lAllen et al.l 
l|2007t l. This latter data set is derived from measurements 
of the baryonic mass fraction in the largest relaxed galaxy 
clusters, which is assumed to be a s tandard quantity for such 
systems (see e.g. lAUen et al.|[2004l . for discussion). 

In General Relativi ty j(a) depends in a no n-trivial way 
on both Qui and w{a) (|Blandford et al.ll2004l ). In general, 
there is no one-to-one mapping between models with con- 
stant j and models with constant w. A powerful feature of 
the standard dynamical approach is that all ACDM models 
have w = —1 which make it easy to search for departures 
from ACDM. Likewise, the use of the jerk formalism im- 
bues the kinematical analysis with a similar important fea- 
ture in that all ACDM models are represented by a single 
value of J = 1. The use of the jerk formalism thus enables 
us to constrain and explore departures from ACDM in the 
kinematical framework in an equally effective manner. More- 
over, by employing both the dynamical and kinematical ap- 
proaches to the analysis of a single data set, we explore a 
wide r set of quest i ons th an w ith a single approac h. We note 
that ISahni et al.] (|2003l ') and I Alam et al.l (|2003l ) also drew 
attention to the importance of the jerk parameter for dis- 
criminating models o f dark ener gy and /or modified gravity. 
IChiba fc Nakamural (|l998l ') and I Caldwell fc Kamionkowskil 
(|2004l ) also showed its relevance for probing the spatial cur- 
vature of the Universe. 

Using the three kinematical data sets mentioned above, 
we find clear evidence for a negative value of qo (current 
acceleration) and a positive cosmic jerk, assuming j con- 
stant. The concordance ACDM model provides a reason- 
able description of the data, using both the new kinemati- 
cal and standard dynamical approaches. We also search for 
more complicated deviations from ACDM, allowing j{a) to 
evolve as the Universe expands, in an analogous manner to 
dynamical studies which allow from time-variation of the 
dark energy equation of state ™(a). Our analysis employs 
a Chebyshev polynomial expansion and a Markov Chain 
Monte Carlo approach to explore parameter spaces. We find 
no evidence for a time- varying jerk. 

This paper is structured as follows; in section [5] we de- 
scribe our new kinematical approach. In section |3] we de- 
scribe the scheme adopted for polynomial expansions of j{a). 
Section 3] includes details of the data analysis. The results 
from the application of our method to the supernovae and 
X-ray cluster data are presented in section [5] Finally, our 
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main conclusions are summarized in section [T] Throughout 
this paper, we assume that the Universe is geometrically flat. 



2 THE KINEMATICAL AND DYNAMICAL 
FRAMEWORKS FOR LATE TIME COSMIC 
ACCELERATION 

2.1 Previous work 

The expansion rate of the Universe can be written in terms 
of the Hubble parameter, H = d/a, where a is the scale 
factor and d is its first derivative with respect to time. The 
current value of the Hubble parameter is the Hubble Con- 
stant, usually written as Hq. Under the action of gravity, and 
for negligible vacuum energy, the expansion of the Universe 
is expected to decelerate at late times. Contrary to this ex- 
pectation, in the late 1990s, type la super novae experiments 
IIRiess et al ] |l998l : |Perlmutter et al.lll999l ) provided the first 
direct evidence for a late time accelerated expansion of the 
Universe. In particular, the present value of the decelera- 
tion parameter, go, measured from the supernova data was 
found to be negative. In detail, the deceleration parameter q 
is defined as the dimensionless second derivative of the scale 
factor 



and in terms of the scale factor, 



q{a) 



H 



{aH)' 



(1) 



(2) 



where the 'dots' and 'primes' denote derivatives with respect 
to cosmic time and scale factor, respectively. 

The current 'concordance' cosmological model, ACDM, 
has been successful in explaining the SNIa results and 
all other precision cosmology measurements to date. To- 
gether with it's theoretical simplicity, this makes the ACDM 
model very attractive. However, as discussed in the intro- 
duction, the concordance model does face significant theo- 
retical challenges and a wide-range of other possible models 
also provide adequat e descriptions of the current data (see 
ICopeland et~aLll2006l . for an extensive review). 

An excellent way to distinguish between models is to 
obtain precise measurements of the time evolution of the 
expansion of the Universe. Given such data, a number of 
different analysis approaches are possible. In searching for 
time evolution in the deceleration parameter, as measured 
by current SNIa data, iRiess et al.l (|2004 ) assumed a linear 
parameterization of q{z), 



dq 

■ go + 

az 



(3) 



These authors measured a change in sign of the deceleration 
parameter, from postive to negative approaching the present 
day, at a redshift zt — 0.46 ± 0.13. Using this parameteri- 
zation for q{z), the definition of the deceleration parameter 
given by equation ((TJ, and integrating over the redshift, we 
obtain that for this model the evolution of the Hubble pa- 
rameter is given in the form 



E{z) = H{z)/Ho = (1 + z)<'+«~''')e'''", 
where g' = dq/dz. 



(4) 



However, since the origin of cosmic acceleration is un- 
known, it is important to recognize that the choice of any 
particular parameterized expansion for q{z) is essentially ar- 
bitrary. Indeed, when (or if) a transition between deceler- 
ated and accelerated phases if inferred to occur can depend 
on the parameterization used. lElgarov fc Multamakil (|2006l ) 
showed that using the linear parameterization described by 
equat ion (|3]) and fitting to the SNIa data set of lAstier et al.l 
(|2005h a transition redshift zt ~ 2.0 is obtained which, un- 
comfortably, lies beyond the range of the data used. 

Transitions between phases of different cosmic accelera- 
tion are more naturally described by models incorporating a 
cosmic 'jerk'. The jerk parameter, j{a), is defined as the di- 
mensionless thi rd derivative of the sca le factor with respect 
to cosmic time (|Blandford et al.ll2004l ) 



./V 1 f'd 



and in terms of the scale factor 



2^2 



(5) 



(6) 



where again the 'dots' and 'primes' denote derivatives with 
respect to cosmic time and scale factor, respectively. 

In such models, a transition from a decelerating phase at 
early times to an accelerating phase at late times occurs for 
all models with go < and a positive cosmic jerk. Note that 
a Taylor expansion of the Hubble par ameter around small 
redshift^s (|Visseij|2004l :[R less et aLll2004 ) contains the present 
value of both the deceleration and jerk parameters, go and 
jo ■ Such Taylor expansions arc inappropriate for fitting high- 
redshift objects (|Blandford et al.„ 2004: Lindc r 200(^), such as 
tho se included in the data sets used here. 

iBlandford elall (|2004l ) describe how the jerk parame- 
terization provides a convenient, alternative method to de- 
scribe models close to ACDM. In this parameterization, flat 
ACDM models have a constant jerk with j{a) — 1 (note 
that this neglects the effects of radiation over the redshift 
range of interest, which is also usually the case when mod- 
eling within the dynamic framework). Thus, any deviation 
from j — 1 measures a departure from ACDM, just as devia- 
tions from w = — 1 do in more standard dynamical analyses. 
Importantly, in comparison to dynamical approaches, how- 
ever, the kinematical approach presented here both explores 
a different set of models and imposes fewer assumptions. The 
dynamical approach has other strengths, however, and can 
be applied to a wider range of data (e.g. CMB and growth 
of structure studies) , making the kinematical and dynamical 
approaches highly complementary. 

It is interesting to note that, in principle, any partic- 
ular dynamical parameter space will have its own phys- 
ical limits. For instance, within the dynamical {Qra,w) 
plane, models with w < —1, known as 'phantom' dark 
energy model s, violat e the dom inant energy condition 
(|Carrofl et a l.| [2003. : Onemh fc Wooda rd 2004) and present 
serious proble ms relatin g to the treatment of dark energy 
perturbations |Hu|[2005| : IVikmanI l2005l : ICaldwefl fc DoranI 
l2005l : IZhao et Z||2005| ) when w{z) crosses the boundary 
w = —1. Current data allow model s with w < — 1 
(IWefler fc Lewis! l2003l: iRiess et al.l l200 4l; I Allen et al.ll2004 
lAstier et al.l l2005l : ISpergel et aH l2006l : ICabre et all |2006| ) 
and models in which w{z) crosses the boundary w = —1 
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Jjassal et al.l2005l: 


Corasaniti et al.ll2004 


Seliak et al.ll2005l 


iRaoetti et alj 2005 


; UDadhve et al.l200a 


Zhao et alJ 2006h 



However, another dynamical parameter space, coming e.g. 
from a different gravity theory, might not pathologically suf- 
fer from such boundaries around the models allowed for cur- 
rent data. 



Since the {qo,j) plane (see below) is purely kinemati- 
cal, i.e. no particular gravity theory is assumed, we are not 
forced to interpret j = 1, or any locus in this plane, as a 
barrier. Note, however, that caution is required in extending 
the results from the kine matical analysis beyond th e range 
of the observed data (see lAmin fc BlandfordI l2007l . for de- 
tails). For example, inappropriately extending a jerk model 
to very high redshifts could imply an unphysical Hubble pa- 
rameter at early times, i.e., these models do not have a Big 
Bang in the past. 



2.3 Standard dynamical framework 

For comparison purposes, we have also carried out a stan- 
dard dynamical analysis of the data in which we employ a 
dark energy model with a constant dark energy equation of 
state, w. From energy conservation of the dark energy fluid 
and the Friedmann equation, we obtain the evolution of the 
Hubble parameter, H{z) — HoE{z), 

E{Z) = [an(l + Zf + (1 - 0^)(1 + ^)3{l + »)]l/2^ (12) 

where Sim is the mean matter density in units of the critical 
density. As with the kinematical analysis, we assume flatness 
and neglect the effects of radiation density. In this frame- 
work, models with a cosmological constant have w = — 1 at 
all times. 



2.2 A new kinematical framework 

For our kinematical analysis, we flrst calculate H{a) given 
j{a;C) where C = (co, ci, cn) is the selected vector of 
parameters use d to describe th e evolu tion of j(a) (see be- 
low). Following iBlandford et"al] (|2004l l we rewrite the defin- 
ing equation for the jerk parameter ([6]) in a more convenient 
form 



a V"(a) - 2j{a)V{a) = 



(7) 



where ' denotes derivative with respect to a and V{a) is 
defined as 



V{a) 



(8) 



We specify the two constants of integration required by (O 
in terms of the present Hubble parameter Ho and the present 
deceleration parameter go as follows 



V{1) 



1 

2' 



V'(l) 



(9) 



where a{to) = 1 at the present time to- Here the first condi- 
tion comes from -ff(l) = Hq and the second from 



-no 



qo ■ 



(10) 



The Hubble parameter, H{a), obtained from equations 
(O, ((8]) and ((9]) is used to calculate the angular diameter 
{cLa) and luminosity (dL) distances for a fiat Friedmann- 
Robertson-Walker-Lemaitre (FRWL) metric 



dA{a) = a d\,(a) = -jr- a \ 
-no J a 



a''E{a) 



■ da. 



(11) 



where c is the speed of light. These theoretical distances 
di^{a) and dA{a) are then used in the data analysis (see sec- 
tion gj. 

Our framework provides a simple and intuitive approach 
for kinematical studies. For models with go < (> 0), 
the Universe is currently accelerating (decelerating) . Models 
with go < and j{a) — 1 (constant) are currently accelerat- 
ing and have the expansion evolving in a manner consistent 
with ACDM. Any significant departure from j = 1 indicates 
that some other mechanism is responsible for the accelera- 
tion. 



3 EVOLVING JERK MODELS 

Our analysis allows for the possibility the cosmic jerk pa- 
rameter, j{a) may evolve with the scale factor. We have re- 
stricted our analysis to the range of a where we have data, 
[ctmin ~ 0.36, Omax = 1]- In Searching for possible evolu- 
tion, our approach is to adopt ACDM as a base model and 
search for progressively more complicated deviations from 
this. We begin by allowing a constant deviation Aj from 
ACDM {j — 1). For this model, it is possible to solve the 
jerk equation (O analytically. Using the initial conditions 
listed in @, we obtain 



Via) 



^/a 
2 



p — u 



a + I — — I a 



2p 



(13) 



where p = (1/2)^(1 + 8j) and u = 2(go + 1/4). Note that 
in the (go,j) plane for 



j < 



qo + 2qi go < -1/4 



(14) 



-1/8 go > -1/4 

there is no Big Bang in the past 0. The models allowed by 
our combined data sets do not cross this boundary. 

For the next most complicated possible deviation from 
ACDM, we have j(a;C) = i**^™ + Ai(a; C). Here = 
1 and j{a; C) is the cosmic jerk for the cosmology in question. 
In order to meaningfully increase the number of parameters 
in the vector C, we employ a framework constructed from 
Chebyshev polynomials. The Chebyshev polynomials form 
a basis set of polynomials that can be used to approximate 
a given function over the interval [—1, 1]. We rescale this 
interval to locate our function Aj{a;C) in the range of scale 
factor where we have data: 



(l/2)(amin + a n 



(l/2)(a„ 



(15) 



where a is the scale factor in the range of interest and Ac 
is Chebyshev variable. The trigonometric expression for a 
Chebyshev polynomial of degree n is given by 



Tn(cic) = cos(n arccos Oc). 



(16) 



These polynomials can also be calculated using the recurrent 
formula 



^ Allowed {qo,j) values are those for which the equation V(a) = 
has no solution in the past (a < 1) jAmin &: Blandfordll2007l ) . 
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Tn+i(ac) = 2acTn(ac) — Tn-i(ac), 



(17) 



where ro(ac) 
are ri(ac) = 



1 and, for example, the next three orders 
, Tiiac) = 2ac' - 1, T3(ac) = 4ac^ 



3ac, 

etc. Using a weighted combination of these components, any 
arbitrary function can be approximately reconstructed. The 
underlying deviation from ACDM can be expressed as 



(18) 



where the weighting coefficients form our vector of pa- 
rameters, C = (co, ci, cn). Thus, using equation (fTS)) 
we produce different parameterizations for increasing A''. 
With higher A'^'s we allow a more precise exploration of the 
[qo, j{a;C)] parameter space. However, it is clear that this 
process will be limited by the ability of the current data 
to distinguish between such models. In order to judge how 
many orders of polynomials to include, we quantify the im- 
provements to the fits obtained from the inclusion of pro- 
gressively higher orders in a variety of ways (see below). In 
general, we find that models with a degree of complexity 
beyond a constant jerk are not required by current data. 

We note that approches other than expanding Aj in 
Chebyshev polynomials are possible, e.g. one could include 
the dimensionless fourth derivative of the scale factor as a 
model parameter. However, since ACDM does not make any 
special prediction for the value of this derivative, we prefer 
to use our general expansion in Aj here. 



is a "nuisance parameter" which contains both the absolute 
magnitude and Hq. 

For the [qo,j{a;C)] parameter space, the luminosity dis- 
tance dL{z;0) is directly obtained integrating the solution 
of the difi^erential equation (O with the definition (|ll|l as 
presented in subsection 12.21 For models using linear param- 
eterization of q(z) and/or dynamical models with Qm and 
w, we plug the equations Q and (|12|) . respectively, into the 
equation describing the luminosity distance for a fiat FRWL 
metric, in units of megaparsecs 



c{l + z) 
Ho 



dz 



E{z-e) 



(19) 



where the speed of light, c, is in km s^^ and the present 
Hubble parameter. Ho, in km(s Mpc)~^. Here the vec- 
tors of parameters for each model are 9 = {qo,dq/dz) 
an d 9 = (Slm i W) re spectively. For the gold sample data 
of iRiess et al] (|2004 l. we use the extinction-corrected dis- 
tance moduli, u°^°(zi) and associ ated errors, af. For the 
SNLS data of (Astier etal] (|2005l ) we use the rest-frame- 
B-band magnitude at maximum light m'^{zi), the stretch 
factor Si and the rest-frame color Ci to obtain //'^"{zi) = 
m'^(zi) — M + a(si — 1) — /3ci. These value s were derived from 
the light curves by lAstier et al.l (|2005h . who also provide 
best-fitting values for a = 1.52 ± 0.14 and /3 = 1.57 ± 0.15. 
For both SNIa data sets, we have 



x'^{9;mo) = ^ 



[/i*(2i;e,Mo)-M°''=(^i;e,M)] = 



(20) 



4 DATA AND ANALYSIS METHODS 
4.1 Type la supernovae data 

For th e analysis o f SNIa data, we use both the 'gold' sam- 
ple of_ ^es^FaL|_j2003) and the first year SNLS sample 
of lAstier et all (|2005D . The former data set contains 157 
SNIa, where a subset of 37 low-redshift objects are in 
common with the data of lAstier et al.l (|2005l ) . lAstier et al.l 
ll2005h contains 11 5 objects. We use the measurements of 
lAstier et all (|2005l 'l for objects in common between the stud- 
ies. Thus, comb i ning b oth data sets we have 120 SNIa from 
the lRiess et al.l (|2004l ) gold sample (157 minus the 37 low- 
redshi ft objects in common) and 115 SNIa from lAstier et all 
(l2005l^ . 

The two SNIa studies use different light-curve fitting 
methods. In order to compare and combine the data, we fit 
the observed distance moduli //'^"{zi) = m°^^{zi) — M, where 
m is the apparent magnitude at maximum light after apply- 
ing galactic extinction, K-correction and light curve width- 
luminosity corrections, and Af is the absolute magnitude, 
with the theoretical predictions, /i"^(zi) = m}^{zi) — M — 
51ogj^Q Di,{zi; 9) + Ho, where _Dl ~ Ho cZl is the i^o-free lu- 
minosity distance, fio = 25 — blogioHo and mo = M + /lo 



^ iRiess et al. I ll2004h presented 16 new Hubble Space Tele- 
scope (HST) SNIa, combined with 170 previously reported SNIa 
from ground-based data. They identified a widely used "high- 
confidence" subset, usually referred to as the gold sample, which 
includes 14 HST SNIa. 

^ 71 SNLS objects, plus 44 previously reported nearby objects. 



where the dispersion associated with each data point, af — 
"^Mi obB ~^ ""hit,! + ""v.i- Here cr^. ^^^^ accounts for flux uncer- 
tainties, crfnt i accounts for intrinsic, systematic dispersion 
in SNIa absolute magnitudes and ; accounts for system- 
atic scatter due to peculiar velocities. The SNLS analysis 
includes an intrinsic dispersion of 0.13104 magnitude^ and 
a peculiar velocity scatter of 300 km/s. The gold sample 
analysis includes 400 km/s peculiar velocity scatter, with an 
additional 2500 km/s added in quadrature for high redshift 
SNIa. 

We marginalise analytically over mo 



f 

J — ( 



r(e) = -21n/ exp 
obtaining 



''{9, mo) ]dmo 



where 



[51ogio Di^{zi;9)-m°'^^izi)f 



(21) 



(22) 



(23) 



. _ V- 51ogio73L(^i;e)-m°fa°(^i) 



SNIa 1 SNIa ' 

Note that the absolute value of = a — {lP/c). For the anal- 
ysis in the s tandard dynamic fram ework, our results agree 
with those of lRiess et ah! (l2004h andlAstier et all (120051'), an d 
the comparison work of Nesseris &i PerivolaropoulosI (|200"5h . 



http: / / snls .in2p3.fr/conf/ papers / cosmol / 
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4.2 X-ray cluster data 

For the an alysis of clu s ter X -ray gas mass fractions, we use 
the data of lAllen et alj (|2007l ). which contains 41 X-ray lumi- 
no us, relaxed g a laxy c lusters, including 26 previously stud- 
ied lAUen et alJ (|2004l ). [Some of the original 26 have since 
been revisited by the Chandra X-ra y observatory lead ing to 
improved constraints (for details see I Allen et alJIioOTl l.] The 
new X-ray data set spans a redshift i nterval 0.06 < z < 1.07. 
Our analysis follows the method of I Allen et al.l (|2004l '). fit- 
ting the apparent redshift evolution of the cluster gas frac- 
tion with the expression 



(25) 



where T = (bQ^Ho^) /[{I + O.wVh) ^m] is the normaliza- 
tion of the /gas(2) curve, dA and (f^^(z) are the angular 
diameter distances {dA ~ dh/iX + •^)^) ^o the clusters for 
a given cosmology and for the reference ACDM cosmology 
(with Ho — 70 km(s Mpc)^^ and fim — 0.3) respectively, 
and Da ~ Ho dA is the _ffo-free angular diameter distance. 
For the kinematical approach we treat the normalization 
as a single 'nuisance' parameter, which we marginalize over 
in the MCMC chains. 

F or the dynamical analysis of the same X-ray data, we 
follow I Allen et all (|2004 ) and employ Gaussian priors on the 
present v alue of the Hubble pa rameter Ho = 72 ± 8 km(s 
Mpc)~^ jFreedman et alj 2001), the mean baryon density 
fibft^ = 0.0214±0.0020 (|Kirkman et aLlbOOSl l and the X-ray 
bias factor b = 0.824zfc 0.089 [determine d from the hydrody- 
namical simulations of lEke et al.l (|l998h , including a 10 per 
cent allowance for systematic uncertainties] . The application 
of these priors leads to an additional constraint on S7m from 
the normalization of the fgas{z) curve. Since the kinemati- 
cal approach does not constrain ilm, the kinematical analy- 
sis does not involve these priors and draws information only 
from the shape of the /gas(z) curve. The dynamical analysis, 
in constrast, extracts information from both the shape and 
normalization. 



4.3 Markov Chain Monte Carlo analysis 

For both the kinematical and dynamical analyses, we sam- 
ple the posterior probability distributions for all parame- 
ter spaces using a Markov Chain Monte Carlo (MCMC) 
technique. This provides a powerful tool for cosmological 
studies, allowing the exploration of large multi-dimensional 
parameter spaces. In detail, we use the Metropolis- 
H astings algorithm imp lemented in the COSMOMC code 
of|L ewis fc Bridle! (|2002l ) for the dynamic formalism, and a 
modified version of this code for the kinematic analysis. 

Our analysis uses four MCMC chains for each combi- 
nation of model and data. We en sure convergence lay ap - 
plying the Gelman- Rubin criterion ijCelman fc Rubinll992h , 
where the convergence is deemed acceptable if the ratio of 
the between-chain and mean-chain variances satisfies R—1 < 
0.1. In general, our chains have R — 1 ^ 0.1. 



http: / /cosmologist.info/cosmomc/ 



4.4 Hypothesis testing in the kinematical analysis: 
how many model parameters are required? 

In the first case, we examined a kinematical model in which 
the deceleration parameter go was included as the only inter- 
esting free parameter [see equation Q with q' = 0]. This is 
hereafter referred to as the model Q. As detailed in section^ 
we next introduced the jerk parameter, j, as an additional 
free parameter, allowing it to take any constant value. We 
refer to this as model J, which has the interesting free pa- 
rameters, go and j. Note that model J includes the set of 
possible ACDM models, which all have constant j — 1. 

We next explored a series of models that allow for 
progressively more complicated deviations from ACDM. In 
each case, the improvement obtained with the introduction 
of additional model parameters, has been gauged from the 
MCMC chains using a variety of statistical tests. In the first 
case, we follow a frequentist approach and use the F-test, 
for which 



F : 



(26) 



where Ax^ is the difi'erence in the minimum x^ between the 
two models, xl is the reduced x^ (x^/^'i where u is number 
of degrees of freedom of the fit, dof) of the final model, and 
Am is the difference in the number of free parameters in the 
two models. Given Am and i/, we calculate the probability 
that the new model would give Ax^ ^ F xt Am by random 
chance. This allows us to quantify the significance of the 
model extension. 

The Bayesian Information Criterion (BIG) provides a 
more stringent model selection cri terion and is a n approxi- 
mation to the Bayesian Evidence (|Schwar j|l978l ). The BIG 
is defined as 



BIC ^ -21nC + k\nN 



(27) 



where £, corresponds to the maximum likelihood obtained 
for a given model (thus, —2 In C is the minimum x^), k is the 
number of free parameters in the model and is the number 
of data points. Values for ABIC < 2 between two models 
are typically considered to represent weak evidence for an 
improvement in the fit. ABIC between 2 and 6 indicates 
'positive evidence' for an improvement, and values greater 
than 6 signify 'strong evidence' for the model with the higher 
BIG (|Jeffrevslll96ll : iKass fc Raftervl Il995l : iMukheriee etal] 
ll998l : lLiddlell2004 ). 

Finally, we have compared the full posterior probability 
distributions for different models, using the Bayes Factor 
to quantify the significance of any improvement in the fit 
obtained. The Bayes Factor is defined Eis t he ratio between 
the B ayesian evidences of the two models (|Kass fc Raftervl 
Il995l l. If P{D\9, M) is the probability of the data D given 
a model M, the Bayesian evidence is defined as the integral 
over the parameter space, 9 



E{M) = P{D\M) 



dep{D\e,M)P{e\M), 



(28) 



where P{9\M) is the prior on the set parameters 9, nor- 
malised to unity. We employ top hat priors for all parame- 
ters and evaluate the integrals using the MCMC samples: 



E{M) 



1 



NA9 



(29) 
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Table 1. The marginalised median values and 68.3 per cent confidence intervals obtained analysing each data set and all three data 
sets together. We show these results for two kinematical models: using only qo (Q model) and extending this parameter space with the 
jerk parameter j { J model). We quote P^r degree of freedom for each model and three different statistical tests to quantify the 
significance of extending the parameter space from Q (go) to J (qoJ)- We quote the difference in Axjq, the probability given by F-test, 
the difference in the Bayesian Information Criterion (BIC) and InBj-Q (where Bj-q is the Bayes factor between the two models). Note 
that combining all three data sots we obtain a significant preference for the J model within all three tests. 



Q model model Improvement 



Data set 


90 


xydof 


90 


j 






F-test 


ABIC 


\nBjQ 


Clusters 


-0.55 ±0.14 


39.6/39 


_Qg-,+0.38 


'^■''^-2.00 


39.6/38 


0.01 


5.6% 


-3.7 


-3.2 


SNLS SNIa 


-0.417 ±0.062 


112.1/113 


-0.65 ±0.23 


1 or, + 1.37 

l.oz_-^ 21 


111.0/112 


1.1 


69.4% 


-3.6 


-2.5 


Gold SNIa 


-0.289 ±0.062 


182.8/155 


-0.86 ±0.21 


o '7c: + l-22 
^•"^-1.10 


174.6/154 


8.2 


99.1% 


3.1 


1.2 


Gold+SNLS+Cl 


-0.391 ±0.045 


300.8/272 


-0.81 ±0.14 




290.1/271 


10.7 


99.8% 


5.1 


3.0 



Table 2. The marginalised median values and 68.3 per cent con- 
fidence intervals obtained analysing all three data sets together. 
We show the results for the constant j model (kinematical) and 
the constant w model (dynamical) and their corresponding 
per degree of freedom. 



Approach Model parameters x^/dof 

Kinematical go = -0.81 ± 0.14, j = 2.16tQ ®5 290.1/271 
Dynamical Un, = 0.306taMi), = -'^■^^to.lt 291.7/272 



where A9 is the volume in the parameter space selected 
to have probability 1 within the top hat priors, N is the 
number of MCMC samples and 9n the sampled parameter 
space. Note that P(D\6n) is the expected probability o f 
the data in the posterior distribution l| Lewis fc Bridle|[200^ . 
The evidence of the model E{M) can be estimated trivially 
from the MCMC samples as the mean likelihood of the sam- 
ples divided by the volume of the prior. It is clear, though, 
that this volume will depend on our selection of the top hat 
priors. In order to be as objective as possible, within the 
Bayesian framework, we use the same priors for parameters 
in common between the two models involved in a compari- 
son. For parameters not in common, we calculate their vol- 
umes subtracting their maximum and minimum values in 
the MCMC samples. 

The Bayes factor between two models Mo and Mi is 
Boi = E(Mo)/E{Mi). If InBoi is positive, Mo is 'preferred' 
over Ml. If InBoi i s negativ e . Mi is preferred over Mo. Fol- 
lowing the scale of IJeflrevsl l|l96lf ). if < InBoi < 1 only 
a "bare mention" of the preference is considered warranted. 
If 1 < In Boi < 2.5, the preference is regarded as of "sub- 
stantial" significance. If 2.5 < ln_Boi < 5 the significance is 
considered to be going from "strong" to "very strong" . 



5 RESULTS 

5.1 Comparison of constant jerk and constant w 
models 

We first examine the statistical improvement obtained in 
moving from the simplest kinematical model Q = [go], in 
which qo is the only interesting free parameter, to model 
J = [go, j(co)], where we include constant jerk j = 1+co (i.e. 
we allow j to take values other than zero). The results ob- 
tained, using the three statistical tests described in subsec- 
tion |4]4] applied to each data set alone and for all three data 
sets together are summarized in Table [1] We find that the 
'gold' sample is the only data set that, on it's own, indicates 
a 'substantial' preference for model J over model Q accord- 
ing to the Bayes factor test. Note that this is not only due to 
the fact that the 'gold' sample extends to higher redshifts, 
thereby providing additional constraining power, but also 
due to the fact that the 'gold' sample hints a small tension in 
the ground-based 'gold' sample data to prefer j > 1 values[f|. 
Combining all three data sets, we obtain a 'strong' prefer- 
ence for model J over model Q, from all three statistical 
tests. Table [T] shows the mean marginalised parameters for 
each model and the la confidence levels. Combining all three 
data sets, we obtain tight constraints on qo = —0.81 ± 0.14 
and j = 2.16lo'7g. result represents the first measure- 
ment of the jerk parameter from cosmological dataQ 

Our dynamical analysis of the same three data sets 
gives w = -1.15t|5 ig and = 0.306l^:^4o (see Table 0. 
Figure [1] shows the constraints for both the kinematical 
{qo,j; top left panel) and dynamical (f2m, w; top right panel) 
models, using all three data sets combined. In both cases, 
the dashed lines indicates the expected range of results for 
ACDM models (i.e. a cosmological constant). We find that 



® An analysis of the 'gold' sample data in which the HST 
supernovae are excluded leads to even stronger preference for 
j > 0: AxjQ = 10.6. In this case, for model J we obtain 
go = -1.17 ± 0.28 and j = i.dbtfll- 

^ Note that Riess et al. (2004) measured jo > at the 2a level, 
where jo comes from a Taylor expan sion of the Hubble parameter 
around small redshifts l lVisseilbOol) . As noted in subsection 12.11 
such an expansion is not appropiate when high redshift data arc 
included, as in the 'gold' sample. 
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Figure 1. A summary of the results from the kincmatical (left panels) and dynamical (right panels) analyses. The top left panel shows 
the 68.3 and 95.4 per cent confidence limits in the (qoj) plane for the kinematical model with a constant jerk, j, obtained using all 
three data sets: both SNIa data sets (Riess et al. 2004; Astier et al. 2005) and the cluster /gas data of Allen et al. (2006). The top right 
panel shows the results in the (Clm^w) plane obtained using the same three data sets and assuming HST, BBNS and b priors. (Note that 
the kinematical analysis does not use the HST, BBNS and b priors). The dashed lines show the expectation for a cosmological constant 
model in both formalisms (j = 1, w = —1, respectively). The bottom panels show the confidence contours in the same planes for the 
individual data sets: the SNLS SNIa data (orange contours), the Riess et al. (2004) 'gold' SNIa sample (blue contours) and the cluster 
/gas data (green contours). Here, the dashed lines again indicate the cosmological constant model. 



both the kinematical and dynamical analyses of the com- 
bined data are consistent with the ACDM model at about 
the 1(7 level. 

It is important to recognise that the results from the 
kinematical and dynamical analyses constrain different sets 
of departures from ACDM. We are using two simple, but 
very different parameterizations based on different underly- 
ing assumptions. The results presented in Figure [1] therefore 
provide interesting new support for the ACDM model. 

The lower panels of figure [1] show the constraints ob- 
tained for the three data sets when analysed individually. 



It is important to note the consistent results from the inde- 
pendent SNIa and X-ray cluster data sets. Note that in the 
dynamical analysis, the X-ray data provide valuable addi- 
tional constraints on Qm, when employing the Ho and Qhh^ 
priors. The overlap of all three data sets in both parame- 
ter spaces highlights the robustness of the measurements. 
Comparing the upper and lower panels of figure [T] we see 
how the combination of data sets significantly tightens the 
constraints. 
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5.2 More complicated kinematical models 

For the combined data set, we have also searched for more 
complicated departures from ACDM by including extra 
model parameters, as described in Section |3l We find no 
significant evidence for models more complicated than a con- 
stant jerk model. In particular, we find a negligible Ax^ be- 
tween models with constant jerk J — [qo,j{co)] and the next 
most sophisticated model Ji — [qo,j{a; co, ci)], and between 
the latter model and the next one, J2 = [qo,j{a; cq, ci, C2)]. 

It is, however, interesting to plot the differences between 
the constraints obtained for each model. Figure [2] shows the 
current la and 2a constraints around the median values 
of j{a) at different scale factors, a, over the range where we 
have data [0.36, 1]. The green, lighter contours show the con- 
straints for the jTi model and the red, darker contours for 
the J model. From this figure it is clear that current data 
provide the best constraints around a ~ 0.77, i.e. z ~ 0.3, 
and that at higher and lower redshift more data are re- 
quired. For the low redshift range, the forthcoming SDSS 
II SNIa data will be helpful. For the high redshift range, 
new HST SNIa and further X-ray cluster data should be 
available in the near future. In the longer term, SNIa data 
from the Large Synaptic Survey Telescope (LSST) [f| and the 
Supernovae Acceleration ProbefSNAP) 0, and X-ray clus- 
ter data from Constellation-X ifj should provide tight con- 
straints on both j{a) and w{a). Future galaxy redshift sur- 
veys covering a high redshift range will also help to tighten 
these constraints, using the baryon osc illation experiment 
(jEisenstein et al.li2005l : ICole et al.ll2005l ). 
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0.4 0.6 0.8 1 

a 

Figure 2. The 68.3 and 95.4 per cent confidence variations about 
the median values for j{a) as a function of the scale factor a, over 
the range where we have data [0.36,1]. Results are shown for the 
constant jerk model (model J7) (red, darker contours) and J7i 
model (green, lighter contours). In both cases, the constraints for 
all three data sets have been combined. The dashed line indicates 
the expectation, j = I (constant) for a cosmological constant 
(ACDM) model. 
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5.3 Comparison of distance measurements 

It is interesting to compare directly the distance curves for 
the kinematical (constant j) and dynamical (constant w) 
models, as determined from the MCMC chains. FigOshows 
the 68.3 and 95.4 per cent confidence limits on the offset 
in distance, as a function of scale factor, relative to a ref- 
erence ACDM cosmology with = 0.27, = 0.73. We 
see that the kinematical and dynamical results occupy very 
similar, though not identical, loci in the distance-scale fac- 
tor plane. For the dynamical analysis, the addition of the 
extra constraint on Qm from the normalization of the /gas 
curve tightens the constraints and pushes the results in a 
direction slightly more consistent with the reference ACDM 
cosmology. 

5.4 Comparison with Riess et al. (2004) 

For comparison purposes, we also present the results ob- 
tained using the linear para meterization o f q(z) described 
by equation ((3]) and used bv I Riess et al.l (|2004l ). Figure |4] 
shows the constraints in the plane {qo,dq/dz) determined 
from each data set, and by combining the three data sets 
(solid, orange contours) . It is clear that the constraints from 
the three independent data sets overlap and that by combin- 
ing them we obtain significantly tighter results than using 
the 'gold' sample alone. 

* http:/ /www. lsst.org/lsst_home.shtml 

^ http://snap.lbl.gov/ 

-"^ http:/ /constellation. gsfc.nasa.gov/ 
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Figure 3. The 68.3 and 95.4 per cent confidence limits on the 
offset in distance as a function of scale factor, relative to the 
reference ACDM cosmology, for both the kinematical (constant 
j; green, shaded curves) and dynamical (constant w; dotted and 
dashed curves) analyses. The dotted curves show the results for 
the dynamical analysis in which the additional constraint on Qm 
from the normalization of the /gas curve is used. The dashed 
curve is for a dynamical analysis where this extra constraint on 
the normalization is ignored. The same MCMC samples used to 
construct Fig [l] have been used. 
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Figure 4. The 68.3 and 95.4 per cent confidence limits in the 
{qo,dq/dz) plane obtained using the SNIa data from the first year 
of the SNLS (Astier et al. 2005) (blue contours), the 'gold' sample 
of Riess et al. (2004) (dashed contours), the cluster /gas data of 
Allen et al. (2006) (green contours) and the combination of all 
three data sets (orange contours). 



6 THE DISTANCE TO THE LAST 
SCATTERING SURFACE 

Finally, we note that there is one further pseudo-distance 
measurement available to us - the distance to the last scat- 
tering surface from CMB data. Although this is not a purely 
kinematical data point, for illustration purposes we show 
the constraints on j{a) that can be achieved if one is will- 
ing to make extra assumptions and include this measure- 
ment. The extra assumptions involved, though strong, are 
well-motivated. In detail, in order to use the distance to 
last scattering, we assume that dark matter behaves like 
standard cold dark matter at all redshifts, an assumption 
well tested by e.g. galaxy cluster, weak lensing and galaxy 
redshift surveys at low redshifts and CMB experiments at 
high-z. We also assume that pre-recombination physics can 
be well described by a standard combination of cold dark 
matter, a photon-baryon fluid and neutrinos, and that any 
early dark energy component has a negligible affect on the 
dynamics. With these assumptions, one can construct the 
comoving angular diameter distance to the last scattering 
surface from cIa = fB(ffldec)/SA, where rs{arec) and Oa are 
the comoving sound horizon at decoupling and the charac- 
teristic angular scale of the acoustic peaks, respectively. For 
a geometrically flat Universe with a negligible early dark 
energy com ponent, we calcula te the sound horizon at de- 
coupling as (|Verde et al.ll2003l ) 

Cs(a) 



3(adGc) 



Ho{0,ma + Sir, 



U/2 



da , 



(30) 



where Cs(a) = c/[l + (3f7ba)/(4f2.^)] is the sound speed in 
the photon-baryon fluid, firad = i^-y + i^v is the present 
radiation energy density, and and Qi, are the present 
photon and neutrino energy densities, respectively. We use 
our X-ray galaxy cluster data, assuming HST, BENS and 



b priors, to determine fl^ = 0.27 ± 0.04 (Allen et al. 2006; 
note that this constraint mainly comes from low-redshift 
clusters). We also use the COBE measurement of the CMB 
temperature To = 2.725 ± 0.002K (|Mather et al.lll999h and 
a standard three neutrino species model with negligible 
masses to obtain firad. For these constraints, we obtain 
rs(2dcc) 146^iOM2a 

From llfinshaw et al.l (|2006l ) we have the multipole of 
the first acoustic peak li = 220.7 ± 0.7. This is related to 
Ia by a shift 0, h = l A(i — </>). Using the fitting formula of 
iDoran fc Liiie^ (|2002h . the BBNS prior for ^ib/ i -^ a s calar 
spectral index = 0.95 ± 0.02 (|Spergel et all |2006| ) and 
assuming no early dark energy, we find 6a = 0.6 ± 0.01 de- 
grees. We then obtain a pseudo-model-independent distance 
to decoupling, d(zdc c) — 13.8 ± l.lGpc, where z^cc = 1088 
(|Spergel et al.ll2006l ). 

Fig [5] shows the tightening of the constraints obtained 
using this "data-point-prior" at high redshift 0. Note that 
figure [S] is plotted on the same scale as figure [5] and shows J 
(red, darker contours) and Ji (green, lighter contours) mod- 
els as before, plus the jTi model (blue contours). Note also 
that here the range of the data is [amin ~ 0.0009, a,nax = 1]. 
Again, using equation (|15|l we rescale the Chebyshev inter- 
val [—1, 1] to locate the functions Aj{a;C) in the range of 
scale factor spanned by the data. The prior information at 
high redshift, from the distance to last scattering, tightens 
the constraints significantly. Evidently, the constraints from 
the kinematic analysis are sensitive to the data quality at 
high redshift. 



7 CONCLUSIONS 

We have developed a new kinematical approach to study the 
expansion o f the history of the Univ erse, building on the ear- 
lier work of lBlandford et al] (|2004l ) . Our technique uses the 
parameter space defined by the current value of the cosmic 
deceleration parameter qo and the jerk parameter j, where 
q and j are the dimensionless second and third derivatives 
of the scale factor with respect to cosmic time. The use of 
this (go, i) parameter space provides a natural framework for 
kinematical studies. In particular, it provides a simple pre- 
scription for searching for departures from ACDM, since the 
complete set of ACDM models are characterized by j = 1 
(constant). 

We have applied our technique to the three best 
available sets of redshift-independent distance measure- 
ments, from type la supernovae studies (|Riess et al.ll2003 : 
lAstier et al.ll2005l ) and measurements of the X-ray gas mass 
fracti on in X-ray l umino us, dynamically relaxed galaxy clus- 
ters (|Allen et al.l |2007| ). Assuming geometric flatness, we 



Note that extending the analysis to the decoupling redshift 
•^dGc = 1088 means that the radiation density becomes non- 
negligible. Although, j can still be calculated as usual, j^CDM 
will not equal 1 at these redshifts. However, the ACDM model can 
then be almo st perfectly described as j^^^'^^ (a) = 1 -|- 2/(1 -|- 
(a/ooq)) (see lAmin fc Blandfordll2007l . for details) where aeq is 
the mean marginalised scale factor at equality, from WMAP data. 
We have explicitly verified that, within the Ic values of Ooq, sys- 
tematic ofi'sets due to the affects of radiation have a negligible 
effect on the derived distances. 
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Figure 5. Icr and 2(t constraints on j{a) over the range (including 
the distance to the last scattering surface) of the data [0.0009,1]. 
Note that this figure and figure [2] are plotted on the same scale 
for comparison purposes. This figure shows the same models as 
figure [2]plus the J72 model, and uses the CMB prior as described 
at the text. The dotted line shows the median j{a) curve for the 
J7i model. 



measure qo = -0.82 ± 0.14 and j = 2.16lo:75 (Figure [T]). 
Note that this represents the first measurement of the cosmic 
jerk parameter, j. A more standard, dynamical analysis of 
the same data gives w — — I.ISIq ^is ^^'^ ~ 0.3061q q4o, 
also assuming flatness and HST, BBNS and b priors (Fig- 
ure [1]). Both sets of results are consistent with the standard 
ACDM paradigm, at about the Icr level. 

In comparison to standard, dynamical approaches, our 
kinematical framework provides a different set of simple 
models and involves fewer assumptions. In particular, kine- 
matical analyses such as that presented here do not assume a 
particular gravity theory. The combination of the kinemati- 
cal and dynamical approaches therefore provides important, 
complementary information for investigating late time cos- 
mic acceleration. The fact that both the kinematical and 
dynamical results presented here are consistent with ACDM 
provides important additional support for that model. The 
fact that the two independent sets of distance measure- 
ments, from X-ray galaxy clusters and supernovae, are indi- 
vidually consistent with ACDM, is reassuring (Figure [T]). 

We have searched for departures from ACDM using a 
new scheme based on the introduction of Chebyshev polyno- 
mials. These orthonormal functions allow us to expand any 
deviation from ACDM, Aj{a;C), as a linear combination of 
polynomials. We use the coefficients of these polynomials, 
C, as fit parameters. The current data provide no evidence 
for a dependence of j on a more complicated than a con- 
stant value. However, higher order terms may be required 
to describe future data sets. In that case, our scheme has the 
advantage that, over a finite interval and using enough high 
order terms, it will provide an acceptable global approxima- 
tion to the true underlying shape. Note that this scheme is 



also applicable to dynamical studies of the evolution of the 
dark energy equation of state, w{a). Note also that Cheby- 
shev polynomial expansions of the same order for w{a) and 
j{a) explore a different set of models. For example, a con- 
stant j ^ 1 model corresponds to an evolving w(a) model 
and vice versa. 

We suggest that future studies should endeavour to use 
both kinematical and dynamical approaches where possible, 
in order to extract the most information from the data. The 
two approaches have different strengths, can be applied to 
with a variety of data sets, and are highly complementary. 
The combination of techniques may be especially helpful in 
distinguishing an origin for cosmic acceleration that lies with 
dark energy (i.e. a new energy component to the Universe) 
from modifications to General Relativity. 
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